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The Quark Gluon String Model (QGSM) reproduces well the global characteristics of the pp 
collisions at RHIC and LHC, e.g., the pseudorapidity and transverse momenta distributions at 
different centralities. The main goal of this work is to employ the Monte Carlo QGSM for description 
of femtoscopic characteristics in pp collisions at RHIC and LHC. The study is concentrated on the 
low multiplicity and multiplicity averaged events, where no collective effects are expected. The 
different procedures for fitting the one-dimensional correlation functions of pions are studied and 
compared with the space-time distributions extracted directly from the model. Particularly, it is 
shown that the double Gaussian fit reveals the contributions coming separately from resonances 
and from directly produced particles. The comparison of model results with the experimental data 
favors decrease of particle formation time with rising collision energy. 
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I. INTRODUCTION 

Experiments at Relativistic Heavy Ion Collider (RHIC) 
have demonstrated that hot and dense matter with par- 
tonic collectivity has been formed in ultrarelativistic 
Au-I-Au colhsions at i/s = 200 ^GeV pj]. Proton-proton 
collisions are conventionally used as a reference to com- 
pare with nuclear collisions and to understand the ob- 
served collective effects. The new interest in general fea- 
tures of pp collisions at ultrarelativistic energies appeared 
after the first publications of Large Hadron Collider 
(LHC) data obtained in pp interactions at y^s = 900 GeV 
and 7TeV 0,13. 

The Bose-Einstein (BE) enhancement in the produc- 
tion of two identical pions at low relative momenta was 
first observed in pp collisions about 50 years ago 
Since then, the developed correlation method, colloqui- 
ally known at present as "femtoscopy technique", was 
successfully applied to the measurement of space-time 
characteristics of the production process at the distances 
of few fermis (1 /m = 10^^^ m) (see, e.g., 043 and ref- 
erences therein). The space-time relative distances are 
"measured" by femtoscopy studies at the points where 
the particles stop to interact. This moment occurs at the 
very late stage of the collision, long after the quark-gluon 
plasma (QGP) or any other exotic state of matter has dis- 
appeared. But signals like the geometric growth of the 
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reaction zone and the specific features of the collective 
flow, generated by QGP pressure gradients, are imprinted 
in the final state as very specific space-momentum corre- 
lations influencing particle spectra and femtoscopic radii. 

The system created in ultrarelativistic pp collisions at 
RHIC and especially at LHC energies can be similar to 
the system created in non-central heavy-ion collisions be- 
cause of the large energy deposited in the overlapping 
region and therefore can also demonstrate collective be- 
havior. The strong argument supporting this point of 
view comes from the observation of the almost identical 
multiplicity and momentum dependencies of the femto- 
scopic radii in pp and Au-|-Au collisions by the STAR 
collaboration at RHIC Q. In particular, the momentum 
dependence of the radii can be linked to the collective 
flow developed in the system . The striking result ob- 
tained by the ALICE collaboration from study of the BE 
correlations in pp collisions at ^/s = 900 GeV [9| is the ab- 
sence of the transverse momentum dependence, whereas 
the event multiplicity dependence holds. 

The aim of the present article is to study hadroniza- 
tion processes in pp collisions at ultrarelativistic ener- 
gies using momentum correlation technique within the 
Monte Carlo quark-gluon string model (QGSM) [l^, ^ 
and to compare results of calculations with the exper- 
imental data obtained at RHIC and LHC. This model 
describes successfully the main characteristics oipp inter- 
actions, such as multiplicity, transverse momentum and 
(pseudo)rapidity distributions in a broad energy range 
from = 200 GeV up to top LHC energy = 7TeV 
fl]| . We try to understand to what extent one is able to 
describe the correlation functions (CFs) in ultrarelativis- 
tic pp collisions within the pure string model picture. 

The paper is organized as follows. A brief description 
of the model features is presented in Sect. HH Special 
attention is given to the concept of the formation time. 
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which plays a very important role for study of the fem- 
toscopy correlations. Section HII] introduces the method 
of correlation functions employed by both the STAR and 
the ALICE collaboration. Model results obtained for pp 
coUisions at ^/s = 200 GeV and ^/s = 900 GeV are pre- 
sented in Sect. II Vl Comparison with the available exper- 
imental data is given as well. The proper choice of the 
baseline is discussed, and ability of the double Gaussian 
fit to restore the contributions of string processes and 
resonances to the correlation functions is demonstrated. 
Finally, conclusions are drawn in Sect. [V] 



II. THE MODEL 

A. Basic features 

Our model is the Monte Carlo realization of the quark- 
gluon string model developed in [T^]. Similarly to the 
dual parton model (DPM) [13], QGSM is based on Gri- 
bov's Reggeon field theory (CRT) 14j] accomplished by 
a string phenomenology of particle production in in- 
elastic hadron-hadron collisions. To describe hadron- 
nucleus and nucleus-nucleus collisions the cascade pro- 
cedure of multiple secondary interactions of hadrons 
is implemented. The model incorporates string frag- 
mentation, formation of resonances, and rcscattcring of 
hadrons. As independent degrees of freedom QGSM in- 
cludes octet and decuplet baryons, octet and nonet vector 
and pseudoscalar mesons, and their antiparticles. The 
momenta and positions of nucleons inside the nuclei are 
generated in accordance with the Fermi momentum dis- 
tribution and the Woods-Saxon density distribution, re- 
spectively. 

Pauli blocking of occupied final states is taken into ac- 
count. Strings in the QGSM can be produced as a result 
of the color exchange mechanism or, like in diffr active 
scattering, due to momentum transfer. The Pomeron, 
which is a pole with an intercept ap{0) > 1 in the 
CRT, corresponds to the cylinder-type diagrams. The 
s-channel discontinuities of the diagrams, representing 
the exchange by n-Pomerons, are related to process of 
2k [k < n) string production. If the contributions of all 
n-Pomeron exchanges to the forward elastic scattering 
amplitude are known, the AGK cutting rules [l^ enable 
one to determine the cross sections for 2fc-strings. The 
hard gluon-gluon scattering and semi-hard processes with 
quark and gluon interactions are also incorporated in the 
model via the so-called hard Pomeron exchange [ll|, [l^ . 
The hard Pomeron is nowadays a standard feature at- 
tributed to a variety of CRT-based microscopic models, 
such as DPM 



[la], PHOJET [17| and EPOS [ig. Its 
presence seems to be necessary to describe the rise of 
multiplicity at midrapidity and pr spectra of secondaries 
in pp interactions at LHC energies within the QGSM [11|] . 
Further details of the MC version of QGSM and its ex- 
tension to A + A collisions can be found in [l^ [111, [l^ . 
Figure [T] displays the pseudorapidity and transverse 
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FIG. 1: (Color online) (a) The charged particle pseudora- 
pidity spectra and (b) their transverse momentum spectra in 
non-single-diffractive events calculated in QGSM for pp colli- 
sions at = 200 GeV (dashed lines), 900 GeV (dash-dotted 
lines) and 7 TeV (solid lines) . Symbols denote the experimen- 
tal data taken from 0, [3, HO] ■ 



momentum distributions of charged particles produced 
in non-single diffractive (NSD) pp collisions at y/s = 
200 GeV, 900 GeV and 7 TeV, respectively. Experimental 
data are plotted onto the QGSM calculations also. Since 
the model reproduces the bulk characteristics of the col- 
lisions quite well, it would be useful to apply the QGSM 
for the analysis of particle interferometry. Note, however, 
that the CRT does not provide the space-time picture of 
the system evolution, thus leaving room for the assump- 
tions concerning the femtoscopy correlations quite open. 
Here one has to rely on approaches developed within the 
framework of the string phenomenology. 



B. QGSM and particle coordinates 

The space-time evolution of the collisions starts from 
the interacting partons, i.e., quarks, diquarks and sea 
quarks distributed randomly in the projectile-target over- 
lapping region. The strings between them are stretch- 
ing and subsequently decaying into hadrons. Due to the 
uncertainty principle it takes time to create a hadron 
from constituent quarks. Also, hadrons are compos- 
ite particles, and this circumstance makes the definition 
of the formation time model dependent. According to 
[23| . two definitions of the formation time or, equiva- 
lently, formation length are eligible in the framework 
of the Lund string model 2J|. In the "yo-yo" case it 
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corresponds to the time/coordinate of the first intersec- 
tion point of the hadron constituents ( "yo-yo" formation 
time) . In the so-caUed constituent case it corresponds to 
the time/coordinate of the point of rupture of the string 
(constituent formation time). In the present version of 
the QGSM the constituent formation time is used. The 
string length L = Mg/ln depends on its mass Ms and on 
the string tension n. The mass of the string is not fixed. 
It is determined by the generation of longitudinal and 
transverse momenta of valence quarks at the string ends, 
that depend on the momenta of colliding hadrons. The 
length of the string varies from the maximum value de- 
termined by the momentum of the incident hadron to the 
minimum value determined by the pion mass. Therefore, 
for the formation of a resonance the mass and length of 
the string must be much larger than for production of a 
pion. 

The formation time t* and coordinate z* of i-th hadron 
in the string center of mass can be expressed via its en- 
ergy E* , its longitudinal momentum p*j and the longitu- 
dinal momenta/energies of all hadrons produced by the 
decay of this string 



(1) 
(2) 



Then we calculate ti in the laboratory frame and 
make the propagation of the coordinates to this point 
{xi,yi,Zi,ti): Qi = aoi + tiPai/E,, a = x,y,z. The ini- 
tial spatial distribution of partons in a proton is found 
to be insignificant for the pion coordinate distributions 
at freeze-out, which are dominated by both the forma- 
tion time of hadrons and decay lengths of resonances. To 
study the possible reduction of the formation time be- 
cause of, e.g., increase of the string tension with rising 
incident energy we introduce in Eqs. (HI)-© the scaling 
parameter a, i.e., k, = an^, where kq = 0.88GeV/fm is 
the default value of the string tension coefficient in the 
QGSM found from comparison with experimental data 
at lower energies jlOi]. The coordinate distributions of 
pions at freeze-out are shown in Fig. [2] for pp collisions 
at -y/i — 900 GeV with a — 1, 1.5 and 3. One can see 
that increase of a makes the coordinate distributions nar- 



III. THE CORRELATION FUNCTION 
REPRESENTATIONS 

The momentum correlations are usually studied with 
the help of correlation functions of two or more par- 
ticles. Particularly, the two-particle correlation func- 
tion CF{pi,p2) = A{pi,p2)/ B{pi,p2) is defined as a ra- 
tio of the two-particle distribution from the same event 
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FIG. 2: The dN/dn, n ^ x {&), y (b), z {c), t (d) distribu- 
tions of pions at freeze-out in pp collisions at ^/s — 900 GeV 
with Q = 1 (dotted line), a = 1.5 (dashed line) and a = 3 
(solid line). 



A{p1tP2) to the reference one. The reference distribu- 
tion is typically constructed by mixing the particles from 
different events of a given class. 

In our simulations the weight of each particle pair is 
calculated according to quantum statistics (QS), using 
particle four-momenta pi and four-coordinates Xi of the 
emission points: w = 1 + cos{q ■ Ax), where q = pi — P2 
and Ax = xi — X2- The CF is defined as a ratio of 
the weighted histogram of the pair kinematic variables to 
the unweighted one. The "ideal" case, CFiiieai{pi,P2) = 
A{pi,p2,'w)/A(jii,P2), uses unweighted pairs from the 
same events as the reference. In experiments one utilizes 
unweighted mixed pairs from different events as the refer- 
ence, namely CFreaUsticiPi,P2) = A{pi,p2,w)/ B{pi,p2). 
There is a difference between the ideal pair distribution 
A(pi,p2) and the mixed one B(pi,p2) due to presence 
of energy-momentum conservation for the pairs from the 
same event and absence of it in pairs from the mixed 
ones. This causes a smooth increase of CFreaiisUc with 
g, which reflects the fact that due to energy-momentum 
conservation the probability of two particle emittence in 
the same direction is smaller than that in opposite di- 
rections. Therefore, a more complex fitting procedure is 
needed for the "realistic CF" than for the "ideal CF". 

Generally, the correlations are measured as a function 
of pair relative momentum four vector q. An invariant 
form of this momentum difference commonly used in the 
one dimensional correlation analysis is (Unv — \/(7q — l^p. 
In both the STAR [8] and the ALICE % experiments the 
correlation function is fitted to a single-Gaussian 

CFs^ngle{q^nv) = [l + A CXp (-i?^^^^g2^^)] D{qinv) , (3) 

where the function D(qinv) takes into account any non- 
femtoscopic correlations including the long-range corre- 
lations due to energy-momentum conservation described 
above. The parameters Rinv and A describe the size of 
pion sources and the correlation strength, respectively. 
Here Rinv is measured in the pair rest frame (FRF). 
Concerning the fit given by Eq.Q we have to note that 
the best way to compare the model simulations with the 
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experimental data is the direct comparison of the corre- 
lation functions. Unfortunately, the CFs are not always 
available and one has to compare the results of the fit, 
that is more complicated. For instance, choice of the 
baseline D(qinv) is rather arbitrary. The baseline should 
describe the CF behavior at large qmv where only the 
conservation laws work, but the region of small qinv re- 
mains terra incognita. Different experiments employ dif- 
ferent extrapolations of the baseline to small qinv, e-g-, 
polynomial extrapolations, EMCIS-fit Monte Carlo 
simulations with PYTHIA and PHOJET 0, that give 
some specific behavior at small qinv due to strong jet con- 
tribution in these models, especially noticeable at large 
kt ■ In order to reproduce the experimental fitting proce- 
dures in a model independent way and make a consistent 
comparison of our simulations with different experiments 
we will use below a flat baseline with D{qinv) — 1 for 
STAR and ALICE data. 

The correlation strength parameter A can differ from 
unity due to the contribution of long-lived resonances, 
particle misidentification and coherence effects. The 
ID correlation functions were studied within the dif- 
ferent ranges of the average pair transverse momentum 
kr ~ Ipt^i +Pt^2|/2 in the mid-rapidity region. 

If large statistics sets are available it is possible to 
perform the 3D correlation analysis. Within realistic 
models, the directional and velocity dependence of the 
correlation function can be used to determine both the 
duration of the emission and the form of the emission 
region, as well as to reveal the details of the produc- 
tion dynamics 0-0] ■ For these purposes the correlation 
functions can be analyzed in terms of the out, side and 
longitudinal components of the relative momentum vec- 
tor q = {qout,qside,qiong} HH, [H . Here qout and q^ide 
denote the transverse components of the vector q, and 
the direction of qout is parallel to the transverse com- 
ponent of the pair three-momentum. The corresponding 
correlation widths are usually parametrized in terms of 
the Gaussian correlation radii Ri 

CF(pi,p2) = 

1 -h Aexp - i?^idc9sidc - ^Lg9io„g) -(4) 

The three dimensional analysis is performed in the lon- 
gitudinal CO- moving system (LCMS), where the pair mo- 
mentum along the beam vanishes. It is possible to com- 
pare the radii measured in LCMS with Rinv by making a 
boost of all radii from LCMS to PRF, namely, RoutPRF — 
IrRout, RsidePRF = Rside, RiongPRF = Riong and aver- 
aging these radn. 

The method used by STAR and ALICE experiments is 
to create a 3D correlation function by filling a three di- 
mensional histogram with the fuU q = {qout,qside, qiong} 
vector in different ranges of the average pair transverse 
momentum kx = |pt^i -1-^^21/2. 





k, = (0.15-0.25) GeV/c 


k, = (0.25-0.35) GeV/c 

1 QGSM,a=1 

n QGSM,a=1.S 

■ QGSM,a=3 




k, = (0.35-0.45) GeV/c 


k, = (0.45-0.6) GeV/c 



' ' o!5 ' 1i„. (Gev;c)) ' ' ' q,„, (Gsv;c) 



FIG. 3: The Tv+n+ CFs for pp at = 200 GeV in four kr 
bins obtained by mixed pair reference distribution. Cuts are 
|77| < 0.5 and 0.12GeV/c j pr i 0.8GeV/c, as in the STAR 
experiment. Calculations are performed with a — 1 (dotted 
line), a = 1.5 (solid line) and a — 3 (dashed line). 
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FIG. 4: The same as Fig|3]but for ^ = 900 GeV/c. Cuts are 
\ri\ < 0.8 and 0.15GeV/c | pr \ l.OGeV/c. Calculations with 
a — 1 (dotted line), a = 1.5 (dashed hue) and a — 3 (solid 
line) are compared to ALICE results 3| with multiphcity M < 
6. 



IV. RESULTS AND DISCUSSION 

The two-pion correlation functions CF(qinv) simulated 
for pp collisions within the QGSM model with the scal- 
ing parameters a = 1, 2 and 3 are shown in differ- 
ent kr ranges in Fig. [3] and Fig. g] for a/s = 200 GeV 
and ^/s = 900 GeV, respectively. The denominator of 
the CF was calculated by means of the mixing proce- 
dure described in Sec. IIIII As expected, smaller forma- 
tion times lead to smaller freeze-out radii of the parti- 
cle sources and, therefore, to larger CFs in the interval 
< qinv < 0.5GeV/c. In Fig. 0] the correlation func- 
tions obtained with the QGSM are directly compared 
to those measured by the ALICE collaboration. The 
ALICE analysis performed for the minimum bias event 
sample gives for the value of the average pseudo-rapidity 
density {dNch/drj) = 3.6, that coincides with the results 
of the QGSM simulations. We compare the QGSM low 
multiplicity sample with the ALICE data at low mul- 
tiplicity bin M < 6. The best description is achieved 
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for the scaling parameter equals to 3. In Fig. |4]one can 
see that the agreement between the shapes of the corre- 
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FIG. 5: The baseline for tt^tt^ correlation functions ex- 
tracted from QGSM calculations (full circles) of pp collisions 
at = 900 GeV in different kr intervals. Solid lines denote 
the fit to polynomial D(qinv) = a + bqinv + cqf„^. 

lation functions calculated within the QGSM and mea- 
sured by the ALICE is rather good till kx < 0.7GeV/c. 
In the last kx bin 0.7 < kr < l.OGeV/c the experimen- 
tal correlation function is about 15% narrower than the 
QGSM one. To understand this effect better the realistic 
correlation functions without QS-weights, i.e. "baselines", 
were constructed in different kx bins as displayed in 
FiglSl The energy-momentum conservation produces the 
long range correlation effects at large qmv , for which the 
calculated values of the CFs lie above the unity. In 
Ref. [Sl the good description of the long range correla- 
tions was obtained within the PYTHIA and PHOJET 
models. In FiglS] the QGSM baseline D{qinv) demon- 
strates complicated behavior qualitatively similar to that 
of the PYTHIA/PHOJET baselines but a bit flatter in 
low qinv interval for large fc^ bins. 
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FIG. 6: One-dimensional ■k'^tt'^ correlation radii as func- 
tions of kr in pp collisions at y/s = 200 GeV (squares) and 
yi = 900 GeV (circles). Open symbols denote STAR [1] and 
ALICE 9] experimental data, full symbols present QGSM 
calculations with a = 1.5 (200 GeV) and a = 3.0 (900 GeV), 
respectively. Both the model results and the data are ob- 
tained from the fit to Eq.Q with the flat baseline. 

Figure [6] presents the fcr dependence of Rmv obtained 
from the fit to Eq. with the flat basehne of the QGSM 



CFs, shown in Fig.[3]and Fig.H The available STAR and 
ALICE data points with flat baseline [1, Q are averaged 
over the multiplicity and compared with the multiplicity 
averaged QGSM correlation functions. The best agree- 
ment with the STAR data 8] was obtained for calcula- 
tions with a — 1.5. 

It was reported in ^ that if PHOJET/PYTHIA base- 
lines are chosen the correlation radii are practically inde- 
pendent on kx within the studied transverse momentum 
range, however, the strength of the kx dependence relies 
heavily on the baseline hypothesis. The ALICE conclu- 
sion about the absence of kx dependence is based on 
the assumption that both PHOJET and PYTHIA cor- 
rectly describe the non-femtoscopic effects at low-gi„„ 
possibly related to minijets. In this case the enhance- 
ment at low-Qiny in the large kx bins is misinterpreted 
as Bose-Einstein enhancement. Since we assign Bose- 
Einstein weights to all pion pairs, we automatically im- 
ply that the enhancement at low-qinv is caused solely by 
the BE correlations. In such a case it will be improper to 
use the PHOJET/PYTHIA or our own QGSM baseline 
to exclude the assumed non-femtoscopy correlations at 
low Qinv ■ The rather successful description of the ALICE 
points within such approach suggests that there are no 
room for non-femtoscopic correlations at low qinv up to 
kx < 0.7GeV/c. 

The ALICE and STAR data points obtained with the 
flat baseline reveal a similar slope in Fig. [51 which is de- 
scribed rather well by the QGSM calculations with the 
scaling factors a = 1.5 and a = 3, respectively. How- 
ever, the higher fc^-bins have larger deviations from the 
experimental points. 
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FIG. 7: (Color online) The space-momentum correlations of 
direct pions produced in QGSM calculated pp collisions at 
\/s = 900 GeV. Line is drawn to guide the eye. 

It is helpful to understand the origin of the strong kx 
dependence of the correlation radii in the QGSM model. 
The LUND hadronization schema described by Eqs. ([1])- 
(Pl) introduces automatically the space-momentum corre- 
lations. The "p-x" correlations for the direct pions dis- 
played in Fig. [7] look similar to the space- momentum cor- 
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relations in hydrodynamic models, where they arise due 
to transverse collective flow. Note that only the parti- 
cles with nearby velocities in their center-of-mass system 
contribute to the correlation function. If the "p-x" corre- 
lations are absent, the whole source is "seen" by the CF 
in any chosen range. Thus, there should be no fcy- 
dependence of the correlation radii. In the presence of the 
"p-x" correlations the particles with close momenta come 
from nearby space regions of the source. Therefore, one 
is measuring not the real geometrical size of the source, 
but rather the size of the regions which emit particles of 
a given momenta, the so-called regions of homogeneity 
[25|. Higher kx pairs should have narrower coordinate 
distributions due to larger "focusing effect". It origi- 
nates from the fact that particles with large momenta fly 
away from each other much quicker than particles with 
small momenta, so in order to be correlated they have 
to be very close in the coordinate space. In Fig. [8lja),(b) 
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FIG. 8: (a) Coordinate distributions of the direct pions in 
PRF in QGSM calculated pp collisions at = 900 GeV 
with a = 3. The three transverse momentum intervals are 
KTl = (0.1-0.25) GeV/c (dotted histogram), KT3 = (0.4- 
0.55) GeV/c (dashed histogram) and KT5 = (0.7-1.0) GeV/c 
(solid histogram). The single Gaussian fit for the KT5 bin 
is shown by the solid line, (b) CFideai for KTl (dotted his- 
togram), KT2 (dashed histogram) and KT3 (solid histogram). 

(c) The same as (a) but for all pions, the single and the dou- 
ble Gaussian fits are shown for KT5 bin by the solid lines. 

(d) The same as (b) but for all pions, double Gaussian fit is 
shown for KT5 bin by the solid line. 

the transverse coordinate distributions are shown in the 
pair-rest-frame together with the corresponding correla- 
tion functions CFideai for the direct pions in three fc^ 
ranges, namely KTl = (0.1 - 0.25)GeV/c; KTi = 
(0.4 - 0.55)GeV/c and KTb = (0.7- 1.0)GeV/c. We 
see that the widths of the Xpup distributions decrease 
with rising kx and the corresponding CFideai become 
narrower. These widths are reproduced within error bars 
by the fit of CFideai to the Gaussian given by Eq. ^ with 

D{qrnv) = 1- 

The important factor influencing the coordinate dis- 
tributions is the ratio of direct pions to pions from reso- 
nance decay. Table |T] presents the fractions of pions from 
decay of the resonances most essentially contributed to 



TABLE I: The fraction of pions from decay of main resonance 
species in QGSM and the path length 1* of these states. 

r (fm) 200 GeV 900 GeV" 



Direct 7r+ 




46.9% 


37.5% 


7r+ from p°'+ — ^ 7r^'*'7r+ 


3.3 


37.1% 


40.7% 


■K^ from u] — > Tr'V^Tr"'" 


28.1 


11.2% 


15.9% 


7r+ from K''^+{K*'0) K-k+ 


8.0 


4.2 % 


5.5 % 



the correlation functions. The path length I* ~ pd/mj^T 
of these states in the cm. frame of two identical pi- 
ons at small value of qiny is listed in Table |T] also. Here 
Pd is the momentum of the decay pion in the resonance 
rest frame [26], m^r is the pion mass and F is the decay 
width. The pions from the decays of rather long-lived 
resonances uj and K* cause appearance of the exponen- 
tial tails in the pion emission function, which distorts the 
Gaussian-like shape of the CF, see Fig. [51[c),(d). Their 
relative contribution decreases with increasing kx due to 
kinematical reasons, whereas the relative contributions 
of direct pions and pions from p decays increase as dis- 
played in Fig. ini This effect also leads to decrease of the 
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FIG. 9: The evolution of the correlation radii Ri„v of identi- 
cal pions with their transverse mass rrit in QGSM calculated 
pp collisions at y^s = 900 GeV. The symbols denote directly 
produced pions (full circles) and pions coming from the de- 
cays of p mesons (full squares), to mesons (open triangles) and 
K* (full triangles), respectively. 

correlation radii with increasing kr- The essentially non- 
Gaussian coordinate distributions that include contribu- 
tions from resonances cannot be fitted well to a single 
Gaussian, however, the double-Gaussian fit reproduces 
its shape properly, see Fig. |Hl[c). By fitting the corre- 
sponding CFs to a single Gaussian one cannot describe 
the narrow peak produced by pions from the resonance 
decays at low qinv ■ On the other hand, using the double- 
Gaussian fitting procedure similar to the one suggested 
in [26,] 

CFdoubieiqinv) = [1 + Ai exp (-i?-„„,ig^™) 
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+ A2 exp {-R^„y^2<lznv)] D{q,nv) , (5) 

where parameters Rinv{i.2) and A(i^2) describe the sizes 
and the correlation strengths of the direct pion source 
and the one of the pions from the resonance decays, re- 
spectively, one gets much better description of the CF 
shape at low qtnv, as shown in Fig. [S^d). 

In order to understand to what extent one is able to 
describe the correlation functions of all particles includ- 
ing the resonances by the different fitting procedures we 
make a comparison of the extracted values of Rmv with 
the Gaussian widths of the coordinate distributions in 
the pair rest frame. The comparison is presented in 
Fig-IHUc) for the "ideal" correlation functions CFideai and 
in Fig. [lU] for the realistic CFs. The extracted parame- 



k, = (0.15-0.25) GeV/c 


^^^^^^k^=^^25-0^4^^ 


^^ooooooooooooooooc 


k, = (0.40-0.55) GeV/c 


k, = (0.55-0.7) GeV/c 






|j^^JO^.O) GeV|c 


o QGSM, a=3 

Double gaussian fit 

Single gaussian fit 
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FIG. 10: The fit of pion correlation functions, obtained in 
QGSM calculated pp collisions at -/s — 900 GeV with a — 
to single Gaussian (dotted line) and double Gaussian (full 
line) in five /ct bins. 

ters are listed in Table |lT] for three ranges, namely, 
KTl = (0.1 - 0.25) GeV/c; KT3 = (0.4 - 0.55) GeV/c 
and KT5 = (0.7 - 1.0)GeV/c. Because of the sharp 
peak of the correlation functions at low qinv the two 
radii restored by the double Gaussian fit vary consid- 
erably. The first one is of the order of 1 fm and has a 
tendency to decrease with rising kx, whereas the second 
one is always larger than 3fm and increases to 13-14 fm 
at high transverse momenta. The second Gaussian is 
quite narrow thus leading to a hair-width difference be- 
tween the single-Gaussian and double-Gaussian curves at 

Qinv 

> O.lGeV/c. 

The ideal 3D correlation functions for ^/s = 200 GeV 
and ^/s = 900 GeV, constructed for the minimum bias 
events and low multiplicity bin, are displayed in Fig. [TT] 
and Fig. [T^l respectively. The calculations were done 
with a = 1.5 and a — 3.0, and the full 3D fit to the 
3D Gaussian given by Eq. ^ was performed. The ex- 
tracted Ri as functions of average kx are presented in 
Fig. [13] and FigfMl One can see that the experimen- 
tal points are rather close to the QGSM ones especially 
for ALICE experimental data, see FiglTil where the low 
multiplicity bin is considered. Note that no integration 
over multiplicity was done in both cases. At 200 GeV 
all radii demonstrate the weak decrease with kx , whereas 
at 900 GeV the radii Rout and Rside are rather flat, the 



TABLE II: Comparison of the Gaussian widths of the co- 
ordinate distributions in PRE shown in Fig[8lc) aXpjip for 
single and double Gaussian fit with the Rinv extracted by us- 
ing different fitting strategies: 1 - "ideal" CE is fitted to the 
single Gaussian Eq.Q with D{qinv) = 1; 

2 - "realistic" CE is fitted to the single Gaussian Eq.® with 
D{qi„v) = 1; 

3 - "realistic" CE is fitted to the single Gaussian Eq.Q with 
D{qinv) = a + bqinv + cq1„^\ 

4 - "realistic" CE is fitted to the double Gaussian Eq.([5]) with 
D{qi„v) = 1; 

5 - "realistic" CE is fitted to the double Gaussian Eq.((5]) with 

D{qi„v) = a + hqinv + cqf„^. 

The selected transverse momentum intervals are 0.1 < /ct < 
0.25GeV/c (KTl), 0.4 < kr < 0.55GeV/c (KT3) and 
0.7 < fcr < l.OGeV/c (KT5), respectively. 



Method 


R 


nvl(2) (fm) 




KTl 


KT3 


KT5 


1 


1.00 


0.77 


0.66 


2 


1.26 


0.84 


0.71 


3 


1.10 


0.84 


0.71 


4 


1.23 


0.81 


0.71 




5.04 


3.26 


13.97 


5 


1.05 


0.81 


0.71 




3.61 


3.25 


13.83 






'^^PRF 


(fm) 


single G. 


3.37 


2.45 


2.96 


double G. 


1.48 


1.08 


1.00 




5.35 


4.72 


4.23 




EIG. 11: Projections of the 3D Cartesian representations of 
the correlation functions onto the qout, qside, and qiong axes, 
for the minimum bias events from pp collisions at 200 GeV for 
four kr ranges. To project onto one g-component, the others 
are integrated over the range < qi < 0.12 GeV/c. 

first point in Rout is lower than the other ones, and only 
Riong demonstrates the decrease with rising kx as was 
observed by the ALIGE collaboration at low multiplicity. 



V. CONCLUSIONS 

The following conclusions can be drawn from our 
study. QGSM calculations show strong dependence of 
the correlation radius on the transverse momentum of a 
pion pair. Similar dependence has been observed by the 
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FIG. 12: The same as Fig llll but for the low multiplicity bin 
Nch < 11 of pp collisions at 900 GeV in five ranges. 
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FIG. 13: Three-dimensional tv'^tt^ correlation radii as func- 
tions of kr in pp collisions at y^s = 200 GeV for minimum 
bias events. Open circles denote STAR experimental data, 
full squares present QGSM calculations with a = 1.5. Both 
the model results and the data are obtained from the fit to 
Eq.lll). 



STAR Collaboration, while the ALICE Collaboration re- 
ported ahxiost constant Rinv with increasing kr- How- 
ever, if the flat baseline is employed instead of the one 
simulated by PYTHIA and PHOJET, the ALICE data 
demonstrate the noticeable fcy dependence as well. The 
origin of such a dependence in the QGSM is traced to the 
space-momentum correlations attributed to microscopic 
string models. If these correlations would be absent, the 
correlation radius Rmv will be independent on the pair 
transverse momentum. 

Pions in the model are produced either directly in the 
processes of string fragmentation or from the decays of 
resonances. The relative contribution of the long-lived 
resonances to pion emission function decreases with ris- 
ing fc^, while the corresponding contributions of direct 



processes and short-lived resonances increase. Therefore, 
the correlation radii of pions also decrease with increasing 
the pair transverse momentum. The fit of the ID corre- 
lation functions to the double Gaussian provides a good 
description of the shape of the CFs at low-qi„„ range and 
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FIG. 14: The same as Fig |13l but for the low multiplicity bin 
iVch < 11 of pp collisions at ^ = 900 GeV. QGSM calcula- 
tions are with a = 3. 



enables us to separate the contributions from the direct 
pions and pions from the resonances. 

It was expected that the size of the freeze-out region 
in pp collisions should increase with rising cm. energy 
from = 200 GeV to = 900 GeV due to the in- 
crease of interaction cross section and the number of 
produced resonances. Surprisingly, the radii measured 
by femtoscopy at y/s — 200 GeV are the same or even 
smaller than the ones at ^/s = 900 GeV, as seen in Fig|6l 
The radii obtained within the standard LUND scenario of 
string breaking and the constituent formation time, im- 
plemented in the QGSM, appear to be larger compared 
to the experimental data. Our analysis favors reduction 
of the formation time with increasing energy of hadronic 
collision. One of the possible solutions is the process of 
string-string interaction via, e.g., fusion of strings that 
leads to increase of the string tension. 
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